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Asymptotic Properties of a Sensor Aiiocation Modei 



D. P. Gaver 
P. A. Jacobs 
M. Youngren 

0. Introduction and Summary 

Sensor resources are to be allocated among 1 non-overlapping but possibly 
contiguous geographical locations called nodes. Each sensor look at a node results 
in an observation, with error, of the number of units that are at the node. In 
military applications units may be individual assets such as tanks or ships, or 
possibly small groups organized as platoons or companies and in geographical 
proximity; in ecological applications, they might be animals, singly or in groups, 
or particular vegetation types. The observations are then to be used to estimate 
the numbers of units at each node. The problem is to allocate sensor resources so 
as to minimize a measure of the error of the estimates, and of the estimate of the 
sum of the numbers at all nodes. It is assumed that the units remain on the nodes 
and do not migrate, although for some sensor types the units must be in motion 
on the node at least some of the time in order for detection to take place. 

The objective of this paper is to propose and study adaptive allocation of sensor 
effort in such a way as to focus the sensor's attention sequentially and 
purposefully on nodes so as to pay most attention to those nodes about which 
the greatest uncertainty, or interest, currently prevails. This problem has features 
in common with adaptive bin-packing problems (cf. Gaver et al., 1995) and 
adaptive allocation of customers to servers (cf. Gaver et al., 1993). 



Numerical examples show that a properly selected sequentially adaptive rule 
will provide estimates of improved precision. 

A related investigation that has recently come to our attention is that by 
Thompson and Seber (1994). Their procedures can also be analyzed using the 
methodology proposed in this paper. 



1. The Model 



There are I nodes. Assume that node i contains r/ units. Sensor resources are 
allocated to one of the 1 nodes at times that occur according to a Poisson process 
with rate A. The sensor resource is allocated to node i with probability aj which is 
tailored to depend on past allocations in a purposeful way. Let Z„(i) be the 
observation of node i; a simple model is that m = 1, 2, ...} are iid with 

binomial distribution with r,- trials and known probability of success pi which is 
the probability of unit detection at node i. We assume for illustration that the 
number of units on a node does not change, although an adaptive scheme of the 
type proposed should effectively follow changes in node population. 

Let Niit) be the number of times in [0,f] that node i is observed by a sensor. 
Let 

V(i;t) = Zi(i) + ... + Zjy.{f)(i) 



be the sum of all the observations during [0,f] for node i. 

Under some conditions V(i; t) has a binomial distribution with Ni(t)ri trials 
and probability of success pi, given Niit), although this assumption may not be 
especially accurate in general. Assuming it to be adequate for the moment, an 
estimate for the number of units at node i based on observations made during 
[0,0 is 




v 0‘;0 

Ni{t)pi 



( 1 . 1 ) 



2 



with 




( 1 . 2 ) 






(1.3) 



There has been much statistical attention paid to estimating the number of trials 
in a binomial distribution, given probability of success; for recent discussion see 
Hall (1994). Here we use the simplest such estimator. 

Now consider the following adaptive allocation rule. A sensor resource arriving 
at time t is allocated to node i with probability 



where the hj are strictly positive sufficiently smooth functions such that 
hjicPx, cVy) = cPhjix, y). 

One possible form for hi is 



if it is assumed that given N,(f), Vi(t) is binomially distributed with mean N,(0pz, 
then from (1.3) it follows that for this function h{, sensor resources tend to be 
allocated to those nodes for which the variance, and hence its square root, the 
standard error of the estimated number of imits, is the largest; the probability of 
allocation to the most uncertain node (by this measure) increases rapidly, 
approaching imity as yinaeases. This tends to bring down that standard error 
quickly, and to equate standard errors of the estimates across nodes. Clearly, 
alternative measures of overall sensor performance are feasible, and possibly 
desirable, such as ones that endeavor to equalize fractional or percent error on 





(1.4) 



'Lhi(Nj(t),Vj{t)) 




(1.5) 
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nodes, or ones that also respond to an independent measure of importance of the 
units on a node. Additionally, node contents may be of various types, which can 
be considered. Thus, the present discussion is of an illustration of an adaptive 
allocation scheme. 

The purposeful allocation of (1.4) introduces dependence between {N,(f); f > 0) 
and {V,(0; t ^ 0}. Asymptotic results for the means E[Ni(t)] and E[Vj(0] as the rate 
of the Poisson process A — > are obtained in Section 2. It is shown that the 
purposeful allocation estimate fj(f) is asymptotically vmbiased. Section 3 presents 
results for the asymptotic means for specific form of function hi (1.5). Section 4 
discusses asymptotic results for the second moments of {Ni(0; t > 0} and {V,<0; 



t > 0} and presents approximate expressions for the Vflr[f,(f)] and Var 



im 



Section 5 describes a simpler Poisson approximation. Section 6 presents results 
from simulation experiments. 



2. First-Moment Calculations and Asymptotics 

Note that for /j > 0 

E[Ni{t + /i)|N(f),y(0] = Ni{t)+Xai{N,V)h. 

Thus 



E[Ni{t + h)] = E[N,(0]+ AE[ai(N,V)/i]. 
Assuming derivatives exist, we have 



Similarly, 

E[Vi{t + /i)|N(0, V(0] = Vi{t) + Vte. 



( 2 . 1 ) 



( 2 . 2 ) 



(2.3) 

(2.4) 
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Thus, 



E[l^(( + h)] = £[^(f)] + ;U»iftE[a,.(N, V)A]. 
Assuming derivatives exist, we have 

|E[V,(t)] = npiAE[a,(N,V)]. 



(2.5) 



(2.6) 



Assume oi is of the form (1.4) where hi is sufficiently smooth with 
hiicPx, cVy) = cVhiix, y) for constants c and hi(x, y)>0 for all x, y > 0. We have 



I'';K'W.Vy(/)) 



(2.7) 



XAy(N/()/A,Vj(()/A)' 

} 

Assume lim N,(f)/A = m,(f), lim V^(f)/A = z;,(t), lim -^E[N,(f)/A] = -^m,(t), 
X^oo x-^oo X-^oodt at 

and jta |E[V;(()/A] = ^n,((). 

Dividing both sides of (2.3) and (2.6) by A and letting X—^oo, the bounded 
convergence theorem yields 



and 



Thus, 



i 



dt ' ‘ ‘ 

i 



Vi{t) = riPimi{t). 



(2.8) 



(2.9) 



(2.10) 
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3. Purposeful Allocation 

The equations for (2.8), can be solved for special functions hj. For 

example, assume 



hi{x,y) = 



Oj ■ y{^-Pi) 

r v2r,2 

/ ^ Pi ] 



for 7>0 andfl,- >0. 



(3.1) 



Since 



Vi{t) = mi{t)riPi, 








In this case (2.8) can be rewritten as 



(3.2) 



where C( = (1 - pi)/pi- 
The functions 

mi{t) = K{oi+riCiY^^^^'^h (3.4) 

with 



‘k+Vi ”i(0 

aj + rjCj THiit) 



K = 




y=i 



(3.5) 



are a solution to (2.8). 

Thus, in the limit as f — > «>, the probability that a sensor resource at time t 
looks at node i is 



«i = 






ri) 



\r/(r+i) ■ 



] 



(3.6) 
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If y=0, then 



ai = lH 

and the allocation is equally likely. 

If oo, then 



(3.7) 



^ - ^+ncj 
j 



(3.8) 



If, further, the probability of detection on node i, p,, is constant for all the nodes, 
then the probabilistic allocation for y— > «<. is roughly proportional to the number 
of units on the node. 



4. Scaling and Approximate Variances of Estimators 



Let 




4x 






(4.1) 

(4.2) 



In Appendix A, moment generating functions are used to show the 
asymptotic normality of {(Xj{t), Vf(f)); f ^ 0) as A — > «•. Further, differential 
equations for their second moments are obtciined. 

Rewriting 

Nj(t) = hnj(t) + SXi(t) (4.3) 

Vy{f) = to^(f) + VIVy(l) (4.4) 

An approximate variance of the estimator of the number of units on node i, 




Mil 

N,(() Pi 



(4.5) 
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can be computed using the "delta method" as follows; cf. Bickel and Doksum 
(1977). To begin. 



Pi 



m 

m): 



A Taylor expansion yields 



v;(0 _ 

Ni(t) mi{t) + Xi{t)l4X 



= Mi + 

m,(f) 



1 

m,(t) 



[m/^] 



Vjit) 

rrtiitf 



(x,(0/Va) 




Thus, 



(4.6) 



(4.7) 




2^E[V;.(0X,.W] 



(4.8) 



= ^^{E[v;^(0] + (wifE[xf(f)]-2riPiE[Xi(f)yj(()]}. 



(4.9) 



Hence, an approximate variance of the estimate of the number of units on 
node I 



Kar[^.(()] 



4T4-rH>i"(0]+(':PifE[x?(0]-2TOE[x,(()v;(0]] 

Pi 2mj(f) ‘ ^ ‘ ^ \ (4,]0) 
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To approximate the variance of the sum of the estimates, a Taylor expansion 
yields 



'm) 


— p 


m) Nj(t)_ 


— I-j 



_ ^ 1 1 









--MM) Efx,(0X,.(()l + c{f 



for I /. 
Thus, 



Coi;(n(f),f;(t)) = —^Cot; 
V } / p.py 



Vi(/) V,(() 



(4.11) 



N,(f)'Ny(() 



(4.12) 



Expressions (4.10) and (4.12) can be used to approximate Var 



Init) 



5. A Simpler Poisson Approximation 

Assume the probability of allocation cq is independent of (V,<f), N,(t)); then 
the number of looks at node i is a Poisson process with rate Aa,- independent of 
the other nodes. Further 



E[Niit)] = l + a^^ 


(5.1) 


E[Vi{t)\ = riPi[l + ait] 


(5.2) 


Var[Ni{t)] = ait 


(5.3) 


Var[Vi{t)] = riPi{l -pi) + a,f[r,p,(l - p,) + 


(5.4) 
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(5.5) 



Cov[Ni{t),Vi{t)] = riPiait 

where we assume that at time 0 each node is looked at once. 
In this case (4.10) becomes 



Var[n{t)\ 



^ 1 m{i- Pi) 
Pi 



(5.6) 



and the estimators fj{t) are independent. 

Assume purposeful allocation is adapted with function hi as in (3.1). A simple 



approximation to Var 



im 



t 



can be obtained by neglecting all covariances and 



assuming the number of looks at node i, {N,(f); t > 0), is a Poisson process with 
rate Aa,- where a/ is determined by (3.6). In this case 



Var 



L i 



“I 



1 rfi-Pi) 



fPi + 



(5.7) 



6. Numerical Examples 

Suppose there are 3 nodes with rj units on node j with r\ = 49, r 2 = 25, and 
T 3 = 16. The probabilities of detecting a unit on node j, py, are pi = 1/11, p2 = 0-5, 
andp 3 = 10/11. 

The variance of the estimate of the sum of the units on all the nodes under 
purposeful allocation with hi as in (3.1) was studied using simulation for y= 0, 1, 
10. Each replication of the simulation begins with one observation at each node. 
The times of arrival of a Poisson process vdth rate 1 are then simulated. A node 
for observation at time t is randomly chosen using probabilities 






N,{() 






ViW-Pi) 

mfpi 



10 



for i = 1, 2, 3 with K{t) the normalizing constant. A binomial observation is 

generated for the node chosen. The simulation has 500 replications. 

Table 1 records the sample mean and square root of the sample variance of 

3 

the sum of the estimates for f = 5, 10, 20, 50 where 

1=1 




m 

m)Pi ' 



Also shown are the sample means of the square root of the sum of estimated 



approximate variances (1.3) 




Vim- Pi) 

mfp} ' 



Table 1 

Estimate of Total Number of Units on All Nodes 
n = 49, f2 = 25, rs = 16; pi = 1/11, p2 = 0.5, pa = 10/11 







Simulation 


Approximation 






i 




i 


Square Root of 
Differential 


Square Root 
of Poisson 


Time 


Y 


Sample 

Mean 


Square Root 
of Sample 
Variance 


Sample 

Mean 


Equation 
Approx Variance 


Approx 

Variance 


5 


0 


89.8 


18.0 




15.0 


13.9 


13.9 




1 


88.1 


15.1 




10.7 


12.3 


10.7 




10 


89.3 


15.2 




10.9 


12.5 


10.4 


10 


0 


90.3 


12.4 




11.9 


10.9 


10.9 




1 


89.8 


9.2 




8.2 


9.1 


8.1 




10 


90.1 


10.2 




8.6 


10.2 


8.0 


20 


0 


90.0 


8.9 




8.6 


8.2 


8.2 




1 


89.5 


7.1 




6.1 


6.7 


6.0 




10 


89.5 


8.3 




6.5 


8.0 


6.1 


50 


0 


89.9 


5.8 




5.5 


5.4 


5.4 




1 


89.9 


3.9 




3.9 


4.2 


3.9 




10 


89.5 


6.2 




4.3 


5.7 


4.2 
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Table 1 also displays results of the approximate variance Var 

L i 

obtained using (4.10) and (4.12) and solving the differential equations (2.8) ant 
the second moment equations in Appendix A, (A.27) - (A.32). The equations 
were solved numerically using the 4 ^h/ 5 th order Runge-Kutta-Fehlberg method 
as implemented in MATLAB; (cf. Math Works, 1992). 



Table 1 also displays the simple Poisson approximation to Var 



im 



of 



(5.7) with the a,- in (3.6). 

Recall that the estimate Var 






is vmbiased for both approximations. The 



true value of is 90. The differential equation approximation is close to the 
square root of the sample variance for all values of /for times 10, 20, and 50. The 
Poisson approximation is close to the square root of the sum of the estimated 
approximate variances; both of these approximations are neglecting the 
covariances induced by the purposeful allocation; these covariances become 
more pronounced as /becomes larger. There is no covariance for /= 0, equally 
likely allocation. Note that the Poisson approximation is conservative. However 
for /= 1, the Poisson approximation is within about 10% of the differential 
equation approximation which incorporates the covariances. However the 
difference is larger for y= 10. 

Note that the square root of the sample variance and the differential equation 
approximation suggest that purposeful allocation with /= 1 yields the smallest 
variance of the estimated sum of the numbers of units on all the nodes. A 
rationale for this suggestion foUows. 

Suppose there are a fixed number of looks K that the sensor can take of all 
nodes and the number of units on each node i, rj, is known along with the 
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probability of detecting a unit on node i, p{. Let be the number of looks the 
sensor gives to node i. If each observation has a binomial distribution 



Var 

j'=l J J=1 






Lagrange multipliers can be used to show that the (approximate) ki, i = 1, 2, 3 

1 

are =[^(l-Pf)/PtT • This solution corresponds to 



that minimize Var 



ho) 

.1=1 



the ccj of (3.6) with a, = 0 and y= 1. Thus, if one is interested in minimizing the 

estimated variance of the sum of the number of units on all the nodes, then one 
should look at node i a number of times proportional to If one 



were interested in minimizing the estimated variance of the estimate of the 
number of units on the node with the greatest number of units then one would 
allocate all looks to that node; this corresponds to the purposeful allocation 
policy of y= 

Tables 2 and 3 present results of the simulation experiment with ri = 49, 
T 2 = 25, = 16 and pj = 0.7, p 2 = 0.8, and ps = 0.9. Table 2 presents the simulation 

and approximation results for the estimate of the sum of units on all the nodes. 
Table 3 presents the simulation and approximation results for the number of 
units on the individual nodes. The differential equation results are close to the 
simulated values for times f = 10, 20, 50. The Poisson approximation also seems 
to be adequate. The Poisson approximation may be doing better in this case 
because the probabilities of unit detection are larger. One source of the 
covariance between the estimators is the possibility that V(i, t) may be 0, in 
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which case node i will not be visited very frequently for y> 0 for the purposeful 
allocation with function hi as in (3.1). 



Table 2 

Estimate of Total Number of Units on All Nodes 
ri = 49, t2 - 25, f3 = 16; pi - 0.7, p2 - 0.8, p^ - 0.9 







Simulation 


Approximation 






i 


JirfB 


Square Root of 
Differential 


Square Root 
of Poisson 


Time 


Y 


Sample 

Mean 


Square Root 
of Sample 
Variance 


Sample 

Mean 


Equation 
Approx Variance 


Approx 

Variance 


5 


0 


89.9 


5.52 


3.70 


3.30 


3.30 




1 


89.4 


6.53 


3.15 


3.11 


3.02 




10 


88.9 


8.55 


3.14 


3.06 


2.98 


10 


0 


90.1 


2.75 


2.84 


2.59 


2.59 




1 


89.9 


2.37 


2.45 


2.46 


2.35 




10 


89.7 


2.46 


2.53 


2.51 


2.35 


20 


0 


90.1 


2.04 


2.07 


1.95 


1.95 




1 


90.0 


1.72 


1.80 


1.90 


1.76 




10 


89.9 


1.87 


1.88 


2.06 


1.80 


50 


0 


90.0 


1.30 


1.32 


1.28 


1.28 




1 


90.0 


1.12 


1.17 


1.29 


1.16 




10 


89.9 


1.26 


1.24 


1.35 


1.21 



14 



50 


20 


1—4 

o 




Time 




H- k 

o 




o 


10 




o 


10 


1-4 


O 


10 




o 




49.0 


49.0 


49.1 


49.0 


49.0 


49.1 


48.9 


49.0 


49.1 


48.5 


48.6 


49.0 


Mean 




Node 1 


o o 


0.88 

(0.85) 


1.08 

am 


1-^ i-^ 

^ ON 


1.32 

(1.30) 


1.68 

(1.66) 


1-^ 1—4 

In In 


1—4 i-a 

^ 00 
^ 1-4 


2.33 

(2.20) 


4.85 

(2.00) 


3.95 

(2.27) 


f5 w 

CX (X 

t— » MD 


Square Root 
Sample 
Variance 
(Diff Eqtn 
Approx) 


0.77 

(0.77) 


0.86 

(0.86) 


1.12 

(1.09) 


(61T) 

61T 


1.32 

(1.33) 


1.75 

(1.66) 


1.61 

(1.63) 


1.81 

(1.81) 


2.39 

(2.20) 


ro ro 

H-4 ^ 


137 

(2.38) 


3.11 

(2.81) 


Sample 

Mean 

(Poisson 

Approx) 


~S} 


24.9 


25.0 


25.0 


24.9 


24.9 


25.0 


24.9 


24.9 


25.0 


24.6 


24.8 


24.9 


Mean 




Node 2 


0.72 

(0.72) 


S O 

ON On 


0.63 

(0.59) _ 


1.08 

(1.08) 


0.97 

(0.96) 


0.98 

(0.90) 


1.50 

(1.44) 


1—4 

ro to 

^ ON 


1.36 

(1.20) 


2.94 

(1.88) 


2.38 

(1.62) 


2.20 

(1.53) 


Square Root 
Sample 
Variance 
(Diff Eqtn 
Approx) 
































0.72 

(0.70) 


0.64 

(0.63) 


(6S0) 

190 


(SOT) I 

on I 


(S6-0) 

860 


0.95 

(0.90) 


1—* 1—4 
^ 


1.32 

(1.25) 


1.30 

(1.20) 


1.89 

(1-70) 


1.69 

(1.59) 


1.68 

(1.53) 


Sample 

Mean 

(Poisson 

Approx) 




16.0 


16.0 


16.0 


15.9 


16.0 


16.0 


15.9 


16.0 


16.0 


15.8 


16.0 


16.0 


Mean 




Node 3 


0.71 

(0.66) 


0.45 

(0.45) 


0.33 

(0.32) 


1.06 

(1.00) 


0.71 

(0.68) 


0.52 

(0.48) 


1—4 

ro to 

S' ^ 


0.96 

(0.89) 


0.71 

(0.64) 


1.97 

(L33) 


1.54 

(LOT) 


1.17 

(0.82) 


Square Root 
Sample 
Variance 
(Diff Eqtn 
Approx) 


0.65 

(0.61) 


0.46 

(0.44) 


0.32 

(0.32) 


0.95 

(0.84) 


0.71 

(0.65) 


0.50 

(0.48) 


1.26 

(1.00) 


0.93 

(0.82) 


0.68 

(0.64) 


l-a 

^ to 


1.13 

(0.99) 


0.87 

(0.82) 


Sample 

Mean 

(Poisson 

Approx) 


— -KJ 
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Table 3 

Simulation: 500 replications 
Estimate of Number of Units on Node 
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APPENDIX A 



In the Appendix we present details of the normal approximation. We follow 
an analytical approach used in a different context in Gaver and Jacobs (1995), and 
in Gaver, Morrison, and Silveira (1993). Let the moment-generating function 
(assumed to exist, otherwise use the characteristic function) be 

Y{0, 0 = E[exp{fflV(f) + )}] 

= E 

Condition on (Niit), Viit)), i e {1, 2, ..., /} to obtain 





' / 




1 

(D 

X 

A 


X9yNy(/ + A) + ^,Vy(( + ft) 

U=> J 





= (1 - ^)exp{m{t) + ^(f)} (A.2) 

i 

where 

= (A.3) 

with Z(0 an observation of the number of units on node i. 

Let /i — » 0 to obtain 

|v'(e,«;() = -'^r(e,4;() + ^iE[“i(N(0-V(())exp{fflV(l) + «V(()}c^'i(^i)} (A.4) 

Scaling 
Let 






(A.5) 



exp 



t/=l }=\ 



(A.1) 
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VAt)-Xvi{t) 



(A,6) 



and let A » 1. 
Let 

then 



«>(ftfc()=£[exp{ex{()+m<))]; 



v(el4X,ll~lX-,t) = <p{0,^) ()exp[ VI[ftn(() + )]} ■ 

Thus, we have the following equation from (A.8) and (A.4) 

= ^(0]} 






(A.7) 

(A.8) 



(A.9) 



I r 

+aXe 

j=l - 



xex 



a,|^w(f) + ^X(f),i;(f) + ^y(f)jexp{fiX’(f)+ ^(f)} 
p[4x[0m{t)+ 

Dividing both sides by exp|VX[ft«(f) + ^(0]} 've obtain 

^(p{e, ^;t) + VI q>{e, t)[0m'{t) + ^'(0] 

= -A<p(0,^;t) 



+aXe 

j=l 






(A.IO) 



% 

+'^-^ai(m{t),v{t))-^Yj{t) + 0 ►exp{flX’(f) + ^(f)} 

j ^ 



IS 



Let 



ai(m(t),v(0) 



/ 



(A.ll) 



where hj(x, y) is a sufficiently smooth function with first order partial derivatives 
and hjl^M’xA^y) = ^hj(x,y). 



Let 







H,(y;0= 






(A.12) 



(A.13) 



H.(x;0= (A.14) 

I i 



Hft(y;0= , (ai5) 

I i 



Note that 
/=l 



(A.16) 




Hi{x;t)Xi{t) + H,(y; f)i;(0 + J,Hij,{x;t)X„{t) + £H,jt(y; 0^(0 

k k 
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Since 



2^a,(Aw,(0 + + = 



(A.17) 



this implies that the summed coefficients of 1/VX, 1/A etc. must individually be 
0 . 

Expression (A. 10) can be rewritten as 



^<p{e, fcf) + VI 9(e, fc()[ftn'(0 + M')] 






I L 



71 7a 2 a 2 a 



X< 






i 






’7a 






7 

-I 

ic=l 






-(f 



where b„(i) = E^Z(i)” j the moment of an observation at node i. 



(A.18) 



Let 



^=0 



(A.19) 



Substituting (A.19) into (A.18) results in the following equation for <pQ 
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0 + 0[^'(0 + ^'(0] 

= -X(po{e,^;t) 






X< 



^•(^•(Q/^f(0) 

'Lhj(mj{t),Vj{t)) 

i 






1 

H — p=- 

71 



Hi(5t;()^w(ft«;()+H;(y;()^«>o(e,«;() 



(A.20) 



/ 

■I 

jt=l 



c/y;t ^hk 



Equating terms of order A^/2^ the terms of order A cancel. The terms of order ^/A 
result in the equation 

cpQ{e,^-t)[em\t)^ fy\t)] 



7 

= I 

1=1 






'Lhj[mj{t),Vj{t)) 

i 



+W,(x;f)^^0o(ft ^;7) + H,(y;f)^^o(^/ 0 



1 

jt=l 






(A.21) 



.Oil 
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The terms of order VX cancel if 



f=iZ^K(0/^/(0) 



(A.22) 



In order for this to occur 



i 



and 



(A.23) 



Thus, 



Z^K(0'^;(0) 

j 



Vi(t) = bj(i)njj(t). 



(A.24) 



(A.25) 



Next look for terms of order 1 in (A.20). 






= S 



fh(mi(t),Vj(t)) 

/ 






H,(a:;() «;') + Hj{y;() A?o(« fc') |»i + 



(A.26) 



I it L 






[e.+iiOte]- 



Equations for the joint moments of {(Xj(t), Yj(t)} can be obtained by 
differentiating (A.26) with respect to {6j) and {<^i) and evaluated at 0=^=0. The 
resulting equations are 
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(A.27) 







2H,(Ar;()£[x,(()^] + 2H,(y;()E[X,(f)Y((()] 



) 

k k 



For l^a 



j^E[Xi(l)Xa{t)] 

= H,(x;()E[X,(()X,(()] + H„(y;()£[Y„(f)X,(l)] 
+H,(x-,t)E[XMXi(t)] + H,(y;f)E[X„(()Y,(f)] 
+lH„t(x;<)E[x*(()X,(/)] + 2:H„t(y;()E[X,(l)yjc(l)] 

k k 

+lHtt(x;()£[X„(j)X*(()] + SHrt(y;()E[X„(l)Yk(()]- 

k k 



£ 

dt 




/ 

+2i,(f)[H^(^;()E[X,(()y,(l)] + H,(y;»)E[v'/(<)]] 



+2i^(<)£H,y(x;f)£[x,.(()Y,(/)]+H,y(y;()£[Y,(()y,{()] 

/ 



(A.28) 



(A.29) 
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For j*k 



|E[n(/)y,(,)] 

= {Hy(x;()£[X;(()y*(()] + H,(y;()£[yy(!)yt{l)]}/’l(/) 
+{Ht(jr;()£[x*(»)yy(0] + Ht(y;f)£[yy(t)V|t(f)]}l’l(£) 
+i^(/)lH,((*;')E[x<(On(<)]+H;<(y;<)E[>»(')v<(()] 

I 

+fci WlHw(^;')£[x<(f)v'/(<)] + H«(y;»)E[n(()v)(f)] 

t 

^£[Xt(f)Vk(0] 

i 

+Hi (a:; f ){£[yt (f )X/c (/)] + (()]y, (£)} 

+H„(y; ()[£[n^(0] + £[Xy(()Vi(')lEl(E)] 

+l{H/i/(>^;0£[xy(()yn(f)]+%(y;<)E[vv(»)Vy(f)]} 

/ 

+X{HK(^;0£[Xy(f)Xt(f)] + %(y;t)£[Xi(()yy(t)])f>l(£) 



(A.30) 



(A.31) 
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For £^k 



^E[X<(0n(')] 

+{Ht(x;()E[X^(()Xt(()] + Ht(y;()E[X,(()n(0]}i^(^) 

+SHf/(^;OE[x,(f)n{()] + H,y(y;()E[y*(()V,(l)] 

/ 

+S{Ht,fef)E[x,(()x<(f)]+%(y;')£[y,(')x,(')]}'^('^) 

i 



(A.32) 
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